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ABSTRACT 

Optical and X-ray observations of the quadruply imaged quasar 1RXS J1131 — 1231 show flux ratio 
anomalies among the images factors of ~2 in the optical and ~3-9 in X-rays. Temporal variability of 
the quasar seems an unlikely explanation for the discrepancies between the X-ray and optical flux ratio 
anomalies. The negative parity of the most affected image and the decreasing trend of the anomalies 
with wavelength suggest microlensing as a possible explanation; this would imply that the source of 
optical radiation in RXS J1131 is ~ l0 4 R g in size for a black hole mass of ~ 10 8 M Q . We also present 
evidence for different X-ray spectral hardness ratios among the four images. 
Subject headings: gravitational lensing — quasars: individual (1RXS J1131— 1231) 



1. INTRODUCTION 

1RXS J1131— 1231 was first cataloged as an X-ray 
source by ROSAT , and was subsequently discovered by 
IS fuse et alJ 1)2003(1 to be a rather spectacular quadruple 
gravitational lens. It is a quasar at a redshift of 0.658, 
lensed by an elliptical (but nearly spherical) galaxy at 
a redshift of 0.295 into four images arranged in what 
ISaha fc Williams! <|2003j) call a long -axis quad. In addi- 
tion, a portion of the host galaxy is lensed into an optical 
Einstein ring. Quadruple quasar lenses such as this are 
useful for lens modelers because of the relatively large 
number of constraints they provide. The models yield 
the mass of the luminous plus dark matter responsible 
for the lensing, as well as more limited information on 
the radial profile of the matter distribution in the lens- 
ing galaxy. In spite of these successes, it has proven dif- 
ficult in many cases to match the flux ratios among the 
images predicted by lens models, to those actually o b- 
served (IMao fc SchneiderllllM IMetcalf fc Zhaol l2002). 

Small-scale structure is often invoked to explain these 
flux r atio anomalies, either in t he form of microlensing by 
stars ((Chang fc Refsdallll979t iSchechter fc Wambsganssl 
12002ft or millilensing by dark matter subhaloes 
with masses of ~ 10 6 M lIMao fc Schneider! IT998I: 
iDalal fc Kochanekl l2002|) . Another potential source of 
flux anomalies is differential extinction in the lensing 
galaxy (where the four beams are furthest from each 
other spatially). However, no significant optical or X- 
ray extinction is expected in the outskirts of an early- 
type galaxy such as the one responsible for the lensing 
in RXS J1131. 

In this letter we report on X-ray observations made 
with the Chandra Observatory at a single epoch, and op- 
tical observations made at Magellan at six epochs over 
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the course of fourteen months. We find that the op- 
tical brightness of the source varied by no more than 
0.3 magnitudes during this time. The optical flux ra- 
tios among the four images are discrepant with the lens 
model by factors of ^2, while the corresponding X-ray 
flux anomalies are factors of ^3-9! In addition, the A 
and B images appear to have significantly different X-ray 
spectral hardness ratios than do the C and D images. We 
briefly discuss some possible explanations for the flux ra- 
tio anomalies and the different spectral hardness ratios. 

2. OBSERVATIONS 
2.1. X-Ray Observations 

RXS J1131 was observed for 10.0 ks on 2004 Apr 12 
(ObsID 4814) with the Advanced CCD Imaging Spec- 
trometer (ACIS) on the Chandra X-ray Observatory. 
Each ACIS chip has 1024 x 1024 pixels and is 8.'3 on a 
side (with a pixel size of 0'.'49) . The point spread function 
(PSF) is both energy-dependent and position-dependent. 
Near the aimpoint, the half-power diameter is about 0'.'8 
at 1 keV, broadening to about 1" at 8 keV. The data 
were taken in timed-exposure mode with an integration 
time of 3.14 s per frame, and the telescope aimpoint was 
on the back-side illuminated S3 chip. The data were 
telemetered to the ground in very faint mode. 

The data were downloaded from the Chandra archive, 
and data reduction was performed using the CIA0 3.2.2 
software provided by the Chandra X-ray Center 7 . The 
data were reprocessed using the CALDB 3.1.0 set of cal- 
ibration files (gain maps, quantum efficiency, quantum 
efficiency uniformity, effective area) including a new bad 
pixel list made with the acis_rvm_tiotpix tool. The 
reprocessing was done without including the pixel ran- 
domization that is added during standard processing; 
this omission slightly improves the point spread function. 
The data were filtered using the standard ASCA grades 
and excluding both bad pixels and software-flagged cos- 
mic ray events. Intervals of background flaring were 
searched for, but none were found. 

The IDL-based s oftware package ACIS Extract v3.79 
<|Broos et a,l .1 12002) was used for subsequent reduction 
and analysis. An image of the X-ray data (see Figure Q 
was constructed by reprojecting the events around RXS 

7 http://asc.harvard.edu 
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J1131 in the 0.5-8 keV energy range using a spatial bin 
size of 0'.'16 8 . Model PSFs were produced for the images 
using the CI AD tool mkpsf at energies of 0.277, 1.4967, 
4.51, 6.4, and 8.6 keV. The 1.4967 keV PSF was used in 
a maximum-likelihood reconstruction image of the data 
(10,000 iterations) in order to determine precise positions 
for each of the four lensed images. 

Small apertures (about 0'/3 in radius) centered on these 
positions were used to extract counts and spectra, and 
the CIA0 tools mkacisrmf and mkarf were used to pro- 
duce response files. ACIS Extract corrected the effective 
area response at each energy based on the fraction of the 
PSF enclosed by the extraction aperture at that energy, 
interpolating from the five model PSFs. The apertures 
enclosed roughly 30% of the PSF at 1.5 keV and roughly 
25% of the PSF at 6.4 keV. These small apertures were 
desirable in order to reduce contamination from the other 
lens images, but image A still suffered some small con- 
tamination from images B and C. To correct for image 
A being in the wings of the PSF of both B and C, five 
extraction regions were placed around B and five around 
C at the same radial distance as A. The averages of each 
set of five regions were used for subtraction from A. The 
contribution of the cosmic X-ray background in the lens 
extraction regions is negligible (roughly 0.005 counts). 

The spectru m of each lens image was fit in XSPEC 12.2 
l|Arnaudlll996ft with a powerlaw model absorbed by two 
components, one fixed at the Galac tic column to RXS 
J113 1 of n H = 3.64 x 10 20 cm" 2 Ipickev fc Lockmanl 
1990) and one allowed to vary. Acceptable fits were 
obtained for each image, with reduced \ 2 °f 0.37, 1.3, 
0.59, and 1.8 for images A, B, C, and D, respectively. 
For each of the images, the additional absorption com- 
ponents were consistent with zero, and the upper ends 
of the 1-fi confidence intervals were 6.0, 1.4, 2.1, and 
4.4 x 10 20 cm~ 2 , indicating similar (and small) absorbing 
columns for each image. The powerlaw indices (with 1-a 
confidence intervals) were 1.21+g;^, 1.23±£j^, 1.63±^, 
and 1.93±gjg, indicating some intrinsic spectral differ- 
ences (i.e., not due to absorption) among the lens im- 
ages. To further quantify these differences, the spectral 
hardness ratio (SR) was defined as the observed photon 
flux (photons cm -2 s _1 ) in the 2-8 keV band to that in 
the 0.5-2 keV band (see Table FJ). To characterize the 
intensity of each image, the unabsorbed power-law flux 
was integrated over 0.5-8 keV. Table fl] lists this flux for 
image B and the flux ratio for the other images. 

2.2. Optical Observations 

The lens was observed at six epochs over the course 
of fourteen months in 2004 and 2005 at the Magellan 
6.5-meter Clay telescope at Las Campanas Observatory. 
The observations made use of the Raymond and Beverly 
Sackler Magellan Instant Camera (MagIC), a 2048 x 2048 
direct imaging instrument with a plate scale of 0'.'069 per 
pixel and a 2.4 arcminute field of view. One epoch of 
observations included imaging in three bands — Sloan 
g' ,r' , and i'; the others were limited to the i' band only. 

8 The satellite continuously dithers in a Lissajous pattern on the 
sky, requiring all images to be reprojected. Standard processing 
produces an image with pixels that are 0'.'492 on a side to match 
the physical CCD pixel size. ACIS Extract produces images with 
pixels matched to the size of those in the model PSFs. 



Chandra Image 




i" 

Magellan Image 
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Fig. 1. — Chandra, Magellan, and model images (top to bot- 
tom) of RXS J1131. The raw Chandra image was convolved with 
a Gaussian with a width ~70% the width of the Chandra PSF in 
order to produce a smoother appearance. The model image po- 
sitions and brightnesses are from the model described in ilij.ll as 
are the predicted source and lens positions, marked with circle and 
diamond-shaped caustic, respectively. Also visible is the critical 
curve. 
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TABLE 1 

X-Ray and Optical Properties of RXS J1131 





X-Ray 


Optical 


Model 


F s a 


1.8 x ict 13 


7.5 x 10~ 14 


1 (+) 


FaIFb 


0.18 ±0.04 


1.10 ±0.155 


1.703 (-) 


F C /F B 


0.27 ± 0.03 


0.47 ± 0.063 


0.962 (+) 


Fd/Fb 


0.06 ±0.01 


0.17 ±0.061 


0.113 (-) 


Fb° 


3.2 x 10 44 


1.3 x 10 44 





SR(A) C 1.05 ±0.28 1.10; 0.86 ±0.13 

SR(B) C 1.14 ±0.09 0.96; 0.94 ±0.13 

SR(C) C 0.72 ±0.10 0.96; 0.82 ±0.13 

SR(D) C 0.58 ±0.16 0.93; 0.77 ±0.13 



a Flux in units of ergs cm 2 s 1 ; corrected for a magnification of 
13.9, as determined from a model of the lens. 

b Luminosity in units of ergs s _1 for a source at z = 0.658 and a 
corresponding luminosity distance of 3840 Mpc. No k-corrections 
have been made. 

c Spectral ratios in the X-ray and optical bands. X-ray ratios are 
defined as the observed photon flux in the 2-8 keV band to that in 
the 0.5-2 keV band. The first of the optical spectral ratios is for 
the Sloan g' band to r' band; the second is for the r 1 to %' bands. 
The optical spectral ratios are given in linear flux units. 



TABLE 2 
Optical Variability of RXS J1131 



Julian Date a Image A b Image B b Image C b Image D b 



2453055 


17.69 


17.70 


18.49 


19.71 


2453058 


17.76 


17.79 


18.69 


19.80 


2453135 


17.56 


17.52 


18.51 


19.95 


2453152 


17.54 


17.58 


18.48 


19.64 


2453376 


17.85 


18.16 


18.82 


19.67 


2453475 


18.04 


18.33 


18.95 


19.73 



a Julian date of the Magellan observation. 

b Sloan AB i' band magnitude, after correction for reference stars 
in the same field of view. 



The seeing on these nights varied from 0'.'4 (on the night 
of multicolor observations) to l'.'O. 

Three images from the night with the best seeing were 
used to produce the pseudocolor optical image shown in 
Figure H The g' band was mapped to blue, while r' was 
mapped to green and i' to red. The color stretch has 
been matched to the square root of the flux to bring out 
the faint Einstein ring. 

The DoPHOT PSF-fitting photometry program was used 
to measure the positions and magnitudes of the four 
quasar components and the lensing galaxy, as well as 
five nearby reference stars. The presence of the Einstein 
ring, which has a red color and thus is especially strong in 
the %' band, has resulted in some additional small uncer- 
tainties, both in astrometric and photometric measure- 
ments, which are difficult to quantify. We estimate the 
astrometric errors at O'.'Ol. 

Th e standard stars 101 207 and RU 149F l|Landoltl 
1992) were used to bring the multi-band observations to 
the AB magnitude system. The transformations from 
Johnson colors to S loan AB colors were taken from 
Fuku gita et alJ l|199fifl . The colors thus obtained for the 
four quasar images are reported in Tabled m linear flux 
units. 



The five field stars were used to calibrate the i' band 
photometry for all six epochs to the same magnitude sys- 
tem as the multi-band observations. After this normal- 
ization, the magnitudes of the 5 stars had rms fluctua- 
tions of between 0.02 and 0.07 magnitudes. The mag- 
nitudes of the four quasar components thus derived are 
presented in Table |21 Because of the effects of the Ein- 
stein ring, the uncertainties are larger for these than for 
the reference stars; we estimate them to be 0.15 mag- 
nitudes. Within this level of uncertainty, the data are 
almost consistent with a steady brightness over time, al- 
though there appears to be a slight dimming trend in the 
last three epochs. 

3. DISCUSSION 
3.1. Modeling the Lens 

Following the lead of iSluse et alJ (|2003T> . we used Kee- 
ton's (2001) Lensmodel software to model the lens as a 
singular isothermal sphere (SIS) with an external shear. 
Including the position of the source as well as the lens po- 
sition and strength and the shear strength and direction, 
there were seven free parameters. We used the Magellan 
positions of the four lens components, which had uncer- 
tainties of O'.'Ol, for a total of eight constraints. We did 
not constrain the position of the lens, in order to allow 
for the possibility that mass may not strictly follow the 
light. The fit yielded a reduced \ 2 01 I- 1 - The lensing 
mass was predicted to be 0'.'14 southeast of the observed 
galaxy position, indicating that our model is not perfect; 
this is typical of such simple lens models. We find the 
lens strength to be 1778, and the shear to be 0.12 in a 
direction 73.3° we st of north. These v alues are similar to 
those reported bv ISluse et aD {2003) of 1782, 0.12, and 
14.8° east of north 9 . 

The magnifications predicted by the best fit model are 
-23.7, 13.9, 13.4, and -1.58 for images A, B, C, and D, 
respectively, with signs indicating image parity. These 
appear as flux ratios in Tabled The Fc/Fb and Fa/Fb 
ratios are low by factors of ^2 in the optical, and by 
factors of ~3-9, respectively, in X-rays. The model rel- 
ative intensities were used to create a simulated image 
as it would appear through the Magellan telescope. This 
image is shown in Figure ^ It is clear that the predicted 
flux ratios are different from those observed. 

3.2. Genuine Optical Anomalies 

The X-ray flux ratios clearly appear to be anoma- 
lous, but one may fairly wonder if another relatively 
simple lens model might fit th e optical data better. 
iKeeton. Gaudi. fc Pettersl l)2003[) use the "cusp relation" 
(which predicts in a model-independent way that the flux 
of image A should be approximately the same as the 
sum of images B and C) to establish convincingly that 
a simple smooth model with an elliptical galaxy cannot 
explain the optical flux ratios in RXS J1131. A highly 
flattened model such as an edge-on disk might explain 
them, but the morphology of the galaxy and the round 
Einstein ring rule out this possibility. 

Our own modeling efforts bore out this conclusion. We 
constrained the fluxes to equal the optical values and 
modeled the lens both as an isothermal ellipsoid with 

9 A 90° offset between the two position angles is due to differing 
sign conventions for the shear term. 
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external shear, and as an isothermal sphere with another 
isothermal sphere off-center, to provide shear. These 
models did not fit nearly as well as our best fit model 
above, with the extra contribution to x 2 coming nearly 
exclusively from the flux constraints. We also tried loos- 
ening the constraints on the positions of the images. This 
did improve the flux fits somewhat, but caused the posi- 
tions to be fit far from their observed values. 

3.3. Quasar Variability 

Another possibility for explaining the differences in the 
X-ray and optical flux ratio anomalies might be to invoke 
temporal variability in the intrinsic output of the quasar, 
since the observations in the two wavebands were made 
at different epochs - though the X-ray observation was 
made about half way between the second and third of the 
optical observations. We have shown directly that the 
optical flux did not undergo any major secular changes 
in the intensity during that year. RXS Jl 131 has a lu- 
minosity which is about midrange (on a log scale) for 
quasars. Therefore it may undergo substantial temporal 
vari ability in both its intensity and spectral slope (see, 
e.g.. iGreen. McHardv. fc Lehtol llQfll . However, a sus- 
tained (> 10 4 s) change in intensity by a factor of ~9 
within a day (the time delay difference between images A 
and B) would be quite unusual (Green et al. 1993). Thus, 
it seems doubtful that temporal variability explains the 
principal flux ratio anomaly in this source. 

3.4. Anomalies Due to Substructure 

As is the case for most anoma lous quadruply lensed 
quasars ijKochanek fe D alai 2003J), the sense of the flux 
anomaly in RXS J1131 is to demagnify the brightest sad- 
dle point image (image A), and possibly to further mag- 
nify the brightest minimum (image B). This discrimi- 
nation by image parity is expected for both micro- and 
millilensing (Schechter & Wambsganss 2002; Kochanek 
et al. 2003). In addition, the fact that the anomaly in 
image A is more severe in X-rays, which originate from a 
smaller region than the optical light, supports microlens- 
ing rather than millilensing. 

To help understand what substructure might do to 
the intensities of the optical and X-ray images, we es- 
timate the ratio of the angular sizes of the emitting re- 
gions near the quasar black hole to the Einstein radius 
of a point object (e.g., a star) in the tensing gala xy (see 
also iMortonson. Schechter. fe Wambsganss 2004) . Ob- 
jects orbiting near the central black hole at radius r sub- 
tend a characteristic angle at the earth of: 



D S 



GM- 



BH 



c 2 D s 



(1) 



where Mbh is the black hole mass, R g is the gravita- 
tional radius of the black hole, and D$ is the angular- 
diameter distance to the source. By comparison, the Ein- 
stein radius of a micro- or millilensing point mass is: 



Be 



UGM mL D LS 
V D L D s c 2 



1/2 



(2) 



where M m L is the mass of the micro- or millilensing ob- 
ject, and Dl and Dls are the lens and lens-to-source 
angular diameter distances, respectively. 



We can define a dimensionless ratio of these quanti- 
ties £, which is related to the degree to which micro- or 
millilensing can occur: 



t-E-Mxur- is. 



( Mbh \ 

Vio 8 AW 



M„ 



M G 



-1/2 



D S D LS 



For RXS Jl 131 the above expression reduces to 



i ~ 3 x 10" 



/ M 



BH 



Rg) V10 8 M Q 



Mr, 



Mr. 



J_ 

0.01 



M 



BH 



Mr, 



('lO ) \ M 



-1/2 



-1/2 



(3) 



(4) 



(5) 



where f3 is the characteristic speed of orbiting or free- 
fall objects around the black hole, and we have taken 
D s ~ 1400, Dl ~ 900, and D LS ~ 865, all in units of 
Mpc. 

Thus, X-ray and optical continuum emission which is 
emitted by the accretion disk within several hundred R g 
of the black hole can be substantially microlensed (see 
eq. (4)). Any broad-line emission features (with j3 ~ 
0.01) could be only marginally microlensed (see eq. (5)). 
In contrast, any narrow- line emission region (with [3 ~ 
10 -3 ) would not be microlensed. In this study, we are 
limited to X-ray and continuum optical emission, both 
of which should be about equally microlensed. 

Therefore, the clear differential in the flux ratio anoma- 
lies between the optical and X-ray bands, factors of ^2 
in the former, and ^3-9 in the latter, presents something 
of a puzzle (see 2J. If, on the other hand, the contin- 
uum optical emission originates farther from the center, 
possibly due to scattering of visible light or reprocess- 
ing of higher energy radiation, then the differential flux 
anomalies between X-ray and optical could be explained 
by microlensing. In this case we can directly estimate 
the size of the optical emission region as ~ 10 4 R g (see 
eq. (4)) for a ~10 8 M black hole. 

Finally, to determine if it is plausible to explain 
a factor of 9 demagnification using microlensing, we 
examined the microlensing sim ulations described by 
iSchechter fc Wambsganss! {2002) . For a saddle-point im- 
age with a magnification of ^20 such as image A, the 
probability of a demagnification a factor of 9 or greater 
ranges from virtually zero for a 100% stellar local pro- 
jected mass density to nearly 17% for a mass density 
made of 10% stars and 90% smooth dark matter. We ex- 
pect that at this distance from the galaxy's center, stars 
would make up about 15-30% of the projected mass, and 
so conclude that it is possible for microlensing to explain 
the X-ray anomaly. 

4. SUMMARY AND CONCLUSIONS 

We have analyzed optical and X-ray images of the quad 
lens 1RXS J1131 — 1231 and find anomalous flux ratios 
among the four images that are different in the opti- 
cal than in the X-ray, with the more extreme anomalies 
being present in the X-ray band. In particular, the ra- 
tio Fa/Fb is a factor of 9.4 ± 1.7 smaller in the X-ray 
band that is predicted from the model image. The ef- 
fects of microlensing in connection with anomalous flux 
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ratios have been discussed extensively in the literature 
(see, e.g., Metcalf & Zhao 2002; Mortonson et al. 2004). 
Above, we discuss why we would nominally expect the 
microlensing of both the continuum optical and the X- 
ray images to be almost the same. 

If the flux ratio anomaly differences between the X- 
ray and optical are ultimately resolved via microlensing, 
then we can turn the argument around and infer the 
approximate dimensions of the optical emission region 
(see, e.g., Mortonson et al. 2004). First, we define / = 
Lx+opt/ L e dd, where / is the fraction of the Eddington 
limiting luminosity that the X-ray plus optical luminosity 
represents. The parameter / incorporates the fact that 
Lx+o P t (~ 5 x 10 44 ergs s _1 ) is less than the bolometric 
luminosity which, in turn, is less than Eddington. The 
mass of the black hole can then be written as {Mbh — 
2 x 10 6 Mo)/ f. Equation (4) can then be recast as: 

^-10% (X) , ( 6 ) 

where r opt is the size of the optical emission region, and 
we have taken £ ~ 3 in order to weaken the microlensing 
significantly (see, e.g., Schechter & Wambsganss 2002; 
Mortonson et al. 2004), and M mL ~ 1 M©. 

Perhaps even more intriguing than the flux ratio 
anomalies is the evidence we have found for different 
spectral hardness ratios in the X-ray band among the 
four images (see Table QJ. To test for the statistical sig- 
nificance of these different hardness ratios, we divided all 
the X-ray events from the region of the Chandra image 
displayed in Figure ^ into two bands, 0.5-2 keV and 2- 
8 keV, and formed a single hardness ratio. This turned 
out to be SR = 0.97 ±0.038. We then computed spectral 
hardness ratios for each of the four individual images (see 
Table nj. Finally, we evaluated the difference of each of 
these spectral hardness ratios from the overall mean and 
divided by the statistical uncertainties in the ratio deter- 
minations. We find that image B has an SR value higher 



than the mean by 1.7a, while images C and D both have 
lower SR values than the mean by 2.4o\ Comparing SR 
values among the individual images, A and B are con- 
sistent with having the same SR, as are C and D. How- 
ever, image B has an SR value which differs from those 
of both C and D by 3.1er. Our more detailed spectral 
analyses of the individual sources (see i)2.1|l shows rather 
convincingly that there is very little absorbing column 
density, and therefore low-energy absorption cannot be 
the cause of the different X-ray spectral hardness ratios. 
Since there are time delays among the four different im- 
ages, one might possibly invoke temporal variability as 
an explanation. While this is possible, the C and D im- 
ages have similar spectral hardness ratios, but a rela- 
tively long time delay of ^80 days, while the B and C 
images have quite different hardness ratios and a rela- 
tively short time delay of only a few hours. So in order 
to explain the spectral ratio anomalies using temporal 
variability, one must invoke a hardening of the source in 
the 80 or so days between D and A, which would last 
at least a day to accommodate the spectral hardness of 
B, followed by a very rapid softening in time for C only 
a few hours later. An extensive monitoring campaign of 
RXS J1131 with the Chandra Observatory would be very 
valuable in terms of better understanding what temporal 
variability actually occurs in this source. 
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